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1. Introduction 

Lattice simulations with light Wilson quarks have been made possible by significant algorih- 
mic advances and increased computer power. For example, the PACS-CS collaboration recently 
reported results for pion masses as small as 156 MeV [1], a value that would have been thought 
impossible to reach a few years ago. 

As much as one appreciates this developement, it also raises the question about the size of 
the lattice spacing artifacts. These become more and more pronounced as one decreases the quark 
mass keeping the lattice spacing fixed. In particular, one expects a sizable impact of the lattice 
artifacts in the regime where the quark mass m is of the order a^A^^^^. In this case the explicit 
breaking of chiral symmetry due to the quark mass and due to the Wilson term are of comparable 
size [2]. 

The appropriate tool to study this regime is Wilson chiral perturbation theory (WChPT) [2, 3]. 
Here we report some of our recently obtained results for pion scattering in WChPT [4, 5]. We 
computed the scattering amplitude and the scattering lengths for all three isospin channels to one 
loop order in the chiral expansion. As a result of the 0{a^) terms the scattering lengths do not 
vanish in the chiral limit, and chiral logarithms proportional to a^lnM^ are present in the one- 
loop results. These contributions can obscure the continuum chiral logarithms and compromise 
the determination of Gasser-Leutwyler coefficients from numerical lattice simulations. Using the 
expression presented here should help in that respect. 



2. Setup 

We consider WChPT for Nf = 2 degenerate quark flavors. In this case the chiral effective 
Lagrangian including the p^,m and term reads [6, 7] 

^LO = ^{d^Ld^L^)-^m{L + L^) + ^C2a^L + L^)\ (2.1) 

where Z contains the pion fields in the usual way and (...) denotes the trace over the flavor indices. 
/, B are the familiar LO low-energy coefficients (LEC) of continuum ChPT, and C2 is an additional 
coefficient associated with non-zero lattice spacing effects. 

We keep the O(a^) correction in and therefore promote it to a LO term in the chiral 
expansion. This is justified and necessary in the regime where m ~ <3^Aq(^j-,, the so-called large 
cut-off effects (LCE) regime. In the generic small quark mass (GSM) regime [8, 9], one assumes 
m ~ aAQ(^j5. Results in this regime are easily obtained by properly expanding the results for the 
LCE regime. 

The sign of C2 determines the phase diagram of the theory [2]. If C2 > there exists an Aoki 
phase where parity and flavor are spontaneously broken [10]. The charged pions are massless in 
this phase due to the spontaneous breaking of the flavor symmetry. Outside this phase the pion 
mass is given by 

Mq = 2Bm - lc2(^ , (2.2) 
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and it vanishes at m = 020^ /B. If C2 < there exists a first order phase transition. All three pions 
are massive for all quark masses, and the pion mass assumes its minimal value 

<,min = 2|c2|fl^ (2.3) 

at m = 0. 



3. Pion scattering at tree-level 

The scattering amplitude is straightforwardly calculated from the Lagrangian (2.1). At tree- 
level we obtain 

A{s,t,u) = ^{s-Ml-2c2a^), (3.1) 

in terms of the three standard Mandelstam variables. 

Setting the lattice spacing to zero we recover the familiar result of continuum ChPT [11]. At 
non-zero lattice spacing the amplitude receives a correction proportional to C2- This modification 
also affects the scattering lengths for definite isospin 7 = and 1 = 2: 

32^(^«"7^''^') ' ''° = "l6^(^° + ^''"')- ^^-^^ 

For a = we recover the well-known results first obtained by Weinberg [12]. For non-zero lattice 
spacings, however, the continuum results are modified in such a way that the scattering lengths no 
longer vanish in the chiral limit. Instead, they assume non-zero values of 0{a^). Consequenty, the 
ratio Qq/Mq is not a constant but has the functional form 

n' A' 

with A JO being a constant and Aqq being of order a^. Hence, the ratio a^/M^ diverges in the chiral 
hmit. This divergence was first predicted by Kawamoto and Smit [13]. However, note that the 
coefficient Aqq is of O(a^) rather than of 0(a). This holds even for unimproved Wilson fermions, 
in contrast to what has been expected earlier [14, 15]. 

The divergence in the chiral limit will only be present if C2 > 0, because only in this case can 
the pion indeed become massless. For C2 < the pion mass cannot be smaller than the minimal 
value in eq. (2.3), resulting in the minimal values 

12 

«0,min = ^5^^l^2|^^^ ao,min = 0- (3-4) 

for the scattering lengths. 

Figure 1 sketches the pion mass dependence of the scattering lengths for the two possible signs 
of C2- It seems feasible that measurements of the scattering lengths will allow determinations of 
C2- Extrapolating the data for to the chiral limit one may directly read off C2 as the value at van- 
ishing pion mass. A practical advantage is that the calculation of Oq does not involve disconnected 
diagrams which introduce large statistical uncertainties. 
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Figure 1: Sketch of the scattering lengths as a function of the pion mass at non-zero lattice spacing. The 
left panel shows l^Ttf^a^ as a function of the tree level pion mass Mq. For C2 < the pion mass cannot 
be smaller than Mq^^^^. Nevertheless, extrapolating to the massless point the scattering length assumes the 
value — 2c2a^, as indicated by the dashed line. For C2 > the pion mass can be taken zero. At this mass 
the scattering length also assumes the value — 2c2fl^, now with the opposite sign. The right panel shows the 
analogous sketch for the 7 = scattering length. 

The 7=1 channel is somewhat special. The scattering amplitude for this isospin channel is 
given by A{t,s,u) — A{u,t,s). The C2 contribution drops out in this difference and the scattering 
length a\ is just as in continuum ChPT [11]. This suggests that the scaling violations in a\ are 
small. 

4. Pion scattering at one-loop 

4.1 Power counting in the LCE regime 

The power counting in WChPT is slightly non-trivial if we take the O(a^) term at LO. 
consists of the terms of 0(p^,m,a^). In order to renormalize the divergencies of the loop diagrams 
one needs higher order counterterms in the chiral Lagrangian. These are, besides the standard 
ones of continuum ChPT [11, 16], the terms of order p^a^,ma^,a^. In addition, terms of order 
p^a,ma,a^ also exist and are formally of lower order. Consequently, these should also be included 
in the chiral expansion. Hence, one-loop calculations in the LCE regime require the following 
terms in the chiral Lagrangian: 

LO : p^, rn, 

NLO: p^a,ma,a^ (4.1) 
NNLO : p"^, p^nijrrP', p^a^,ma^,a'^ 

The standard NLO terms of continuum ChPT appear here at NNLO, a consequence of the fact that 
m ~ a^Ap^p in the LCE regime. 
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Figure 2: One-loop diagramms contributing to the four-point function. The diagrams with one or two 
vertices stemming from the O(a^) term in the chiral Lagrangian will give rise to chiral logs proportional 
C2a^MQ\nMQ and (c2fl^)^lnMo. 



4.2 The pion mass to one loop 

The modifications due to extra chiral logarithms are already seen in the one-loop result for the 
pion mass. In terms of the tree-level mass Mq we find [7] 



Ml = Ml 



1 M2 fMl\ 5 Icia^^ fMl\ , 
IH ^— ?^ln H ^^^In --^]+ki 



■■3 

IcjWoa^ , {lc2a^Y 

The coefficients A3, S3 and ^1,^3,^4 are (combinations of) unknown LECs. Eq. (4.2) reproduces 
the well-known result in the naive continuum limit, but is modified for non-zero a, in particular 
by a chiral logarithm proportional to C2a'. Moreover, the coefficient in front of the chiral log is 
much larger than the coefficient in front of the continuum chiral log. Consequently, already small 
values of cjo^ may dilute the continuum log and the characteristic curvature one usually looks for 
in numerical lattice data might be significantly diminished. 

We remark that the ETM collaboration has measured a negative value for C2 in their twisted 
mass simulations. The data [17] for a 0.086 fm and M.„;± k, 300 MeV results in —2c2a^ f« (185 
MeV)2, determined from the pion mass splitting. Such a value completely suppresses the chiral 
log for pion masses around 400 MeV, a value not unusual in present day lattice simulations. 

4.3 The 1 = 2 scattering length at one loop 

The four diagrams in figure 2 contribute to the four-point function at one loop. The vertices in 
these diagrams can be either a vertex also present in continuum ChPT or the vertex proportional to 
C2a^ stemming from the 0{a^) term in Eq. (2.1). The latter give rise to additional chiral logarithms 
cici'MqIuMq and {c2a^)^lnMl in the scattering amplitude and the scattering lengths. 

As an example we give the one-loop result for the 1 = 2 scattering length (for the other isospin 
channels see Ref. [4]): 

,2_ Ml f Ml fl Ml 4 .,\ 2c2a' f Ml 

2c2a^ f 2c2fl' f 1 1 1 \ 



7^. 



-h r2 ■ (4-2) 



167r/2 V 16712/2 [2 jU' 

The parameters K2j comprise all the analytic terms through 0(a2). Therefore, these ai^e constants 
for fixed a and serve as fitting parameters in the analysis of numerical data. As a short hand notation 
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we introduced 

lfn=h+lh-\h + \ (4.4) 

for the combination of LECs entering the scattering lengths a^} 

As expected, (4.3) reproduces the continuum result for a ^ 0. For non-zero a, however, it 
differs significantly from the coiTcsponding one in continuum ChPT. The scattering length does not 
vanish in the chiral limit, but rather seems to diverge as o^XnM^. However, if C2 < no divergence 
is present since the pion mass is bounded by the minimal value in eq. (2.3). If C2 > massless 
pions are in principle possible, but for small pion masses the chiral expansion eventually breaks 
down when C2a^ InM^ becomes of order unity. In that case higher order terms leading to powers 
{c2a^lnM^)" ,n = 2,3,... become relevant and a summation of all these terms is necessary.^ This 
should be kept in mind if one keeps lowering the quark mass at a fixed lattice spacing: Eventually 
the lattice spacing corrections will be enhanced by the unphysical divergence o'XnM'^. 

Result (4.3) contains five unknown parameters: the continuum parameters / and l^^^l as well 
as K21, K"22,C2, i.e. three more than the continuum result. Consequently, the pion mass dependence 
of the scattering length on the lattice can be very different in contrast to what one may expect from 
continuum ChPT. Attempts to fit lattice data using expressions from continuum ChPT may easily 
fail. 

The ratio a^/Ml has the functional form 

an Ann 9 TO 9 ~ InM^ 

^ = ^+Aio+A2oMl+A3oMl\nMl+A4olnMl+A4o^. (4.5) 

The first two terms on the right hand side correspond to the tree level result in eq. (3.3). The 
constants Aqq —A40 represent the five independent fit parameters, while ^40 is not independent. 
The first three terms in (4.5) have already been used in analyzing numerical lattice data [18], but 
the data could not be fitted well. It might be interesting to repeat the analysis with the full result 
in (4.5), even though the data was obtained for heavy pion masses between 500 MeV and 1.1 GeV, 
and ChPT is not expected to be applicable. 



5. Conclusions 

For pion scattering with Wilson fermions we essentially find two modifications due to a non- 
zero lattice spacing. First, the 7 = 0,2 scattering lengths do not vanish in the chiral limit, and 
second, additional chiral logarithms proportional to appear in the one-loop results for these 
quantities. Ignoring these modifications and using the results of continuum ChPT in the analysis of 
numerical lattice data is potentially dangerous, depending on the size of these extra contributions. 
In particular, the extra contributions proportional to the lattice spacing become more and more 
relevant the smaller the quark mass is. This should be kept in mind when one pushes to the physical 
point at non-zero lattice spacing. 

' We emphasize that (4.3) is expressed in terms of the decay constant in the chiral limit, and not in terms of /^r. If 
the latter had been used the LEC I4 [11] would also appear in (4.4). 

^This resummation can presumably be done along the lines in Ref. [7], where a resummed pion mass formula has 
been derived. 
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We conclude with a comment on the coefficient C2, which determines the phase diagram of the 
theory. The 7 = 2 scattering length may provide a handle to obtain an estimate of C2. A practical 
advantage is that the calculation of Aq does not involve disconnected diagrams which introduce 
large statistical uncertainties. We expect that at least the sign of C2 is easily accessible, and this is 
what matters for the phase diagram. 
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